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We study the equilibration of isospin degree of freedom in intermediate energy 
heavy-ion coUisions using an isospin-dependent BUU model. It is found that 
there exists a transition from the isospin equilibration at low energies to non- 
equilibration at high energies as the beam energy varies across the Fermi en- 
ergy in central, asymmetric heavy-ion collisions. At beam energies around 55 
MeV/nucleon, the composite system in thermal equilibrium but isospin non- 
equilibrium breaks up into two primary hot residues with N/Z ratios closely 
related to those of the target and projectile respectively. The decay of these 
forward-backward moving residues results in the strong isospin asymmetry in 
space and the dependence of the isotopic composition of fragments on the N/Z 
ratios of the target and projectile. These features are in good agreement with 
those found recently in experiments at NSCL/MSU and TAMU, implications of 
these findings are discussed. 
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With the availabihty of high intensity radioactive beams at several existing and proposed 
facihties, the isospin degree of freedom in nuclear reactions can be studied in large domains 
of beam energies and projectile-target combinations. In recent experimental studies of the 
isotopic composition of intermediate mass fragments from central collisions of ^°C/, ^^Ar 
and ^°Ca with ^^Fe and ^^Ni at Eheam/A=25, 35, 45 and 53 MeV at NSCL/MSU and TAMU 
m-H, it was found that at Ebeam/A=25 and 35 MeV the isotopic ratios ^Bef^Be, ^^B/^°B 
and ^^C/^'^C increase linearly as a function of the increasing {N/Z)cs ratio of the combined 
target and projectile system, but are independent of the N/Z ratio of the target or projectile. 
This observation strongly indicates the establishment of isospin equilibration in composite 
systems formed in these reactions before the emission of fragments. In fact, previous studies 
on similar systems at a much lower energy showed that the isospin degree of freedom was 
one of the fastest to equilibrate . The most striking and unexpected feature was observed 
from the isotopic ratios in central collisions at Ebeam/A = 53 MeV. It was found that the 
isotopic ratios depend on the N/Z of the target and projectile for entrance target-projectile 
combinations having the same {N/Z)cs ratio, such as '^^Ca+^^Fe and '^^Ar+^^Ni. Moreover, 
very forward and backward angles also show isotope ratios that are not simply a 
function of the {N/Z)cs- It was found that light fragments at backward angles have a much 
greater dependence on {N/Z)target, while at forward angles have a much greater dependence 
on (N/Z) projectile, and could not be explained simply by the pre-equilibrium emission. These 
features clearly demonstrate that the isospin degree of freedom is not equilibrated in the 
reactions at i^beam/^ = 53 MeV on the time scale of emission of the fragments. 

The above observation has some deep implications on the reaction mechanism leading 
to multifragmentation. It not only establishes a relative time scale of multifragmentation 
in these reactions but also indicates that the assumption of isospin equilibrium taken for 
granted in various statistical models for nuclear multifragmentation at intermediate energies 
must be improved at least. To our best knowledge, the prominent features observed in these 
experiments were neither predicted nor observed previously in any model simulations. On the 
contrary, some models which are able to explain many other aspects of multifragmentation 
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and reaction dynamics of intermediate energy heavy-ion collisions are unable to explain these 
features. In particular, statistical models based on the assumption of isospin and thermal 
equilibrium of the composite system fail to show any entrance channel effect as one would 
expect. Calculations using an intranuclear cascade code ISABEL ||^ show that the N/Z of 
the residue is very close to the N/Z of the initial combined system and thus also failed 
to reproduce the features observed at Ebeam/A = 53 MeV. Although the exact origin of 
this failure is not so clear to us, one expects that the reaction dynamics at these energies 
should go beyond the nucleon-nucleon cascade and invoke the nuclear mean field. Therefore, 
the experimentally observed transition from isospin equilibration to non-equilibration as the 
beam energy increases from below to above the Fermi energy remains an interesting question 
for theoretical studies. 

In this Letter we report on results of a study on the equilibration of isospin degree of 
freedom using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (BUU) transport model 
for heavy-ion collisions at intermediate energies The isospin dependence comes into the 
model [^] through both the elementary nucleon-nucleon cross sections ai2 and the nuclear 
mean field U. Here we use the experimental nucleon-nucleon cross sections with the explicit 
isospin dependence 0. The isospin dependence resides in the fact that the cross section 
of neutron-proton collisions is about three times that of neutron-neutron or proton-proton 
collisions at energies interested here. 

The nuclear mean field U including the isospin symmetry term is parameterized as 

Uip, r,) = a{p/po) + bip/poT + (1 - r.)Vc + C^^^^r.. (1) 

Po 

Here po is the normal nuclear matter density, p, pn and pp are the nucleon, neutron and 
proton densities, respectively, equals 1 or -1 for neutrons or protons, respectively, and 
Vc is the Coulomb potential. Corresponding to other, probably more complete forms of 
Skyrme forces, other forms of parameterization for the mean field are possible |IO|-p!2[]. Here 
we use the so called soft nuclear equation of state with the nuclear compressibility K = 200 
MeV. The symmetry term comes from averaging over the constituent two-body forces with 
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Heisenberg components proportional to (ii ■ tj) |T^. Its strength C can be deduced from 
experiments (e.g. nuclear symmetry energies, optical potentials for nucleon scatterings, 
excitations of analogue states in (p,n) reactions). However, the strength deduced varies 



significantly from reaction to reaction and also depends on the energy of the nucleon p3 
The parameter C used in transport model calculations also scatters within a large range 
TT| , P^JT5| . In the present study, we use a constant C = 32 MeV. 



First, to study the energy dependence of the reaction mechanism we show in Fig. 1 the 
density contours in the reaction plane at t= 200 fm/c in head-on collisions of Ar+Ni at 
Ebeam/A = 25, 35, 45 and 55 MeV. The solid contours with p = po/8 essentially bound 
the composite systems or heavy residues formed in the reactions. The dotted contours with 
p = 0.05po are the representatives of free nucleons mainly from pre-thermal-equilibrium 
emissions. The most interesting feature in Fig. 1 is the formation of two heavy residues in 
the reaction at Eheam/A = 55 MeV. While by following the evolution of the heavy residues 
until 300 fm/c, we found no breakup happening to the residues formed in the reactions at 
Ebeam/A = 25, 35 and 45 MeV. 

Thermal equilibrium of the heavy residues can be examined by studying the quadrupole 
moment Qzz defined as 

QUt) = J |^(2P' -Pl-Pl)f{r,m- (2) 

Where f{r,p,t) is the Wigner function from the BUU model calculations. It is clear that 
Qzz = is a necessary but not a sufficient condition for thermal equilibrium. We found 
that the quadrupole moments of heavy residues formed in the reactions at Ebeam/A = 25, 
35 and 45 MeV have been damped to about zero by the time of about 300 fm/c. Since 
the heavy residue formed in the reaction at Eheam/^ = 55 has already broken up into two 
pieces at about 200 fm/c, it is interesting to examine more closely the evolution of Qzz in 
this reaction. This is done in Fig. 2 where Qzz of the heavy residue scaled by its mass 
number A^es is shown as a function of time. For comparison we also show in this figure the 
evolution of Qzz for the reaction at Ebeam/A = 45. It is seen that the residues are very close 



to thermal equilibrium at about 200 fm/c indicated by the very small, decreasing amplitude 
of the oscillation around zero. In the case of Eheam/^ = 55 MeV, the oscillation is shifted 
towards a small positive value indicating a slight longitudinal collectivity or transparency. 
This finding is actually in very good agreement with the transparency found recently at 
Ganil for the same reaction system at similar beam energies using the INDRA detector JIB 



Are the heavy residues observed above in isospin equilibrium ? To answer this question 
we show in Fig. 3 the numbers of neutrons and protons in the residues on the left {Z < 0) 
and right {Z > 0) side of the origin in the head-on reactions of Ar+Ni at Eheam/A =25, 35, 
45 and 55 MeV. The solid lines are the numbers of protons from the projectile, while the 
dotdashed lines are the numbers of neutrons from the projectile which comes in originally 
from the left. The dashed lines are the numbers of protons from the target, while the dotted 
lines are the numbers of neutrons from the target which comes in originally from the right in 
the center of mass frame. It is seen that the numbers of neutrons and protons on both sides 
undergo damped oscillations. The damping is mainly due to nucleon-nucleon collisions and 
particle emissions. While the oscillations are due to both the restoring force of the mean 
field and nucleon-nucleon collisions. At E'beam/^ = 25 MeV the neutron and proton numbers 
on the two sides have become very close to each other and the amplitude of oscillation is 
becoming rather small by the time of 300 fm/c. It indicates that the heavy residue is very 
close to the isospin equilibrium characterized by the space-time independent proton and 
neutron number distributions. At E'^eamM = 35 MeV the oscillation is damped at a higher 
rate to a distribution of particles close to the isospin equilibrium. The isotopic composition 
of fragments emitted from the residues of these low energy reactions after about 300 fm/c 
would therefore essentially reflect the {N/Z)cs of the initial composite system, and there is 
little forward-backward asymmetry. These features are in good agreement with those found 
in the preliminary data of Eheam/^ =25 and 35 MeV 

At higher energies, such as Ebeam/^ = 45 and 55 MeV, there is little oscillation across the 
overlapping region between the target and projectile. This is mainly because the incoming 
momenta of projectile-nucleons and target-nucleons are too high for the mean field and 



nucleon-nucleon collisions to be effective enough during a relatively short reaction time to 
reverse the directions of motion of many nucleons. As a result, there exists a large isospin 
asymmetry or non-equilibration at these two energies. In particular, on the left side of the 
origin the N/Z ratio of the residue is more affected by that of the target while on the right 
side it is more affected by that of the projectile. However, it should be stressed that the N/Z 
ratios on the two sides are not simply those of the target and projectile but a combination 
of the two depending on the complicated reaction dynamics. In the case of Eheam/^ — 
55 MeV, at the time of about 200 fm/c the heavy residue has broken up into two pieces 
with some longitudinal collectivity. The forward moving residue has an average N/Z ratio 
of about 7.2/6.0 and an excitation energy of about 8.6 MeV/nucleon. While the backward 
moving residue has an average N/Z ratio of about 21/18 and an excitation energy of about 
6.8 MeV/nucleon. 

The proper decay of these primary hot residues is beyond the scope of the BUU model. 
Nevertheless, one can imagine that the decay of these forward-backward moving hot residues 
with the different N/Z ratios formed in reactions at higher energies will result in a strong 
isospin asymmetry in space. It is also quite understandable that the isotopic composition 
of fragments emitted from these residues will have much greater dependence on the N/Z 
ratios of the target and projectile, but not much on the {N/Z)cs of the initial composite 
systems. Indeed, this expectation is further evidenced in our systematic studies on the 
system ^^S -\-^^ Zn, ^^Ca-\-^^ Fe and ^°Ti -|-^^ Cr at the same beam energies. These features 
are in good agreement with the indication of the experimental data at £^5eam/^ = 53 MeV. 
Furthermore, the properties of the primary residues calculated in the above are useful as 
inputs for a hybrid model study to compare with the data more quantitatively. Although 
there is no global isospin equilibration in reactions at higher energies, it may not be such a 
bad approximation to assume the establishment of isospin equilibrium in the two primary 
residues formed at E'^eam/^ = 55 MeV. One can thus couple the BUU model with statistical 
models to compare more quantitatively with the experimental data. These studies are now 
in progress and results will be published elsewhere. 
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In summary, we have studied the question of isospin equihbration in intermediate energy 
heavy-ion coUisions using the isospin-dependent BUU model. We found that there exists a 
transition from the isospin equihbration to non-equihbration as the beam energy increases 
from below to above the Fermi energy in central, asymmetric heavy-ion collisions. At beam 
energies around 55 MeV/nucleon, the composite system in isospin non-equilibrium breaks 
up into two forward-backward moving hot residues with N/Z ratios more closely related to 
those of the target and projectile on the time scale of thermal equilibrium. These features 
are in good agreement with that found in recent experiments at NSCL/MSU and TAMU. 
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FIGURE CAPTIONS 

Fig. 1 Density contours in the reaction plane at t=200 fm/c for head-on coUisions of 
Ar+Ni at E^beamM —25, 35, 45 and 55 MeV, respectively. The sohd contours are for 
p — Pq/S, while the dotted hues are for p = 0.05po- 

Fig. 2 The evolution of the quadrupole moment per nucleon of heavy residues formed in 
the head-on collisions of Ar+Ni at E^beamM — ^5 and 55 MeV. 

Fig. 3 The numbers of neutrons and protons in the residue with p > pq/8 on the left 
{Z < 0) and right {Z > 0) side of the origin in the head-on reactions of Ar+Ni at 
Ebeam/A =25, 35, 45 and 55 MeV respectively. The solid hues are the numbers of 
protons from the projectile, while the dotdashed lines are the numbers of neutrons 
from the projectile. The dashed lines are the numbers of protons from the target, 
while the dotted lines are the numbers of neutrons from the target. 
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http://arXiv.org/ps/nucl-th/9505015vl 



This figure "figl-2.png" is available in "png" format from: 
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